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ABSTRACT 

The self averaging properties of the conductance g are explored in Random Resistor 



Networks (RRN) with a broad distribution of bond strengths P(g) ~ g^ ■ The RRN 



O 

in 
o 

■^J- ' problem is cast in terms of simple combinations of random variables on hierarchical lat- 



tices. Distributions of equivalent conductances are estimated numerically on hierarchical 
lattices as a function of size L and the distribution tail strength parameter fi. For net- 



works above the percolation threshold, convergence to a Gaussian basin is always the case, 

O " 

except in the limit fx — > 0. A disorder length is identified, beyond which the system 
is effectively homogeneous. This length scale diverges as £d ~ iy is the regular 

percolation correlation length exponent) when the microscopic distribution of conductors 
is exponentially wide (ix — > 0). This implies that exactly the same critical behavior can be 
induced by geometrical disorder and by strong bond disorder with the bond occupation 
probability p «-> fx. We find that only lattices at the percolation threshold have renor- 
malized probability distributions in a Levy-like basin. At the percolation threshold the 
disorder length diverges at a critical tail strength fx c as \fx — fi c \~ z with z ~ 3.2 ±0.1, a 
new exponent. Critical path analysis is used in a generalized form to give the macroscopic 
conductance in the case of lattices above p c . 
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Scaling properties of transport quantities in random media have attracted great in- 
terest for some time. Electric conduction on random networks, in particular, is of great 
relevance because it may be regarded as a simple qualitative archetype of other, more 
complex and less well understood problems in disordered systems. Such physical systems 
include the classical problem of flow in porous mediafl], diffusion in a random environ- 
ment [2] and even the quantum problem of hopping conduction in dirty semiconductors [3] 
or metal insulator composites [4], among others. These systems can be modeled on a ran- 
dom network where individual bonds represent a 'resistor' whose properties depend on the 
quantity which flows through it. In the case of flow in porous media this resistor represents 
a microscopic permeability or permeance (in analogy with conductance to indicate that it 
may not be an intensive quantity). For hopping conduction, on the other hand, the resistor 
involves the computation of the overlap integral between impurities, and quantum effects 
are taken into account [3]. The physics involved in the computation of the resistors deter- 
mines an associated probability distribution function (PDF) for values on the equivalent 
resistor network. The resulting distribution for hopping conduction is exponentially wide 
(r = roexp(-e) with e, say, uniformly distributed), while in the case of flow through porous 
media in the capillary regime it has been argued that the distribution has a power law 
component and a large permeability cutoff [5] . In the same fashion, continuum percolation 
results in a random network with a power law distribution of resistance P(r) ~ l/r M+1 [6]. 
In addition to the statistical disorder discussed above, one can also introduce geometri- 
cal disorder by allowing voids in the network with a certain probability. Networks with 
a broad distribution of bond values at the percolation threshold have been extensively 
studied numerically and by field-theoretical techniques [7] [8] [9]. 

Once the distribution of the elementary resistors is set, as discussed above, one is 
interested in computing the macroscopic properties of such networks. How does the prob- 
ability distribution of the equivalent conductance 'renormalize' as the network is rescaled? 
What is the effective support for transport, given the elementary distribution of resistors? 
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In brief how do the random variables combine and thus 'interact' with the microscopic 
probability distribution (MPD) to give a macroscopic result. In this work we explore the 
scaling properties of conduction on hierarchical networks where the problem can be sim- 
ply and explicitly stated in terms of a non-linear combination of random variables. This 
approach is very attractive since it dwells on the problem of stable limit distributions for 
combinations of random variables (other than additive) for which there are few rigorous 
results. This approach has been adopted recently by Derrida and Griffiths[10], Halpin- 
Healy[ll] and Roux et al[12] in the context of directed polymers in a random medium 
and by Bouchaud, Le Doussal and Georges in a series of papers reviewed in ref.[2] in the 
context of anomalous diffusion. 

The paper is organized as follows: In section I we discuss hierarchical lattices and how 
resistor network composition laws are easily implemented. We then show the results for the 
evolution of the PDF as a function of system size for hierarchical networks at, and above the 
percolation threshold. The existence of two basins of attraction is obtained as a function of 
the MPD tail strength parameter. In section II we discuss the existence of a disorder length 
£d below which strong heterogeneities persist, and derive the conductance at this scale by 
critical path analysis introduced by Ambegaokar, Halperin and Langer (AHL) [3] [13]. This 
length scale determines the crossover distance over which the system self averages. We also 
show the divergence of £d as one approaches P(g) ~ l/g and determine the associated 
critical exponent numerically and using simple analytical arguments. We find that this 
exponent is exactly that of regular percolation, and establish a mapping between the 
probability p of regular percolation and the tail strength parameter \x. In section III we 
point out the parallelism between Levy limit theorems and limit distributions found in 
random resistor networks. Analogies to other disordered systems are also discussed. 

I. Hierarchical Lattices. 

Hierarchical lattices are self similar structures generated by an iterative procedure that 
produces the lattice at order m + 1 from order m, by substituting every bond by a chosen 
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motif (see fig.(l)). With this procedure one obtains lattices of different effective dimensions 
d e and connectivities, which have proven useful in exposing the qualitative behavior of low 
dimensional systems[14]. Fig.(l) shows the so called 'Berker lattice' of effective dimension 
d e = 2 and the Arcangelis, Redner and Coniglio (ARC) lattice[15] frequently used as a 
model for percolating backbones (d e = log 4/ log 3.) 

Bonds on the hierarchical lattices are assigned conductances g, independently chosen 
from the probability distribution 



The ranges assigned to the conductance g secure that the distribution is normalized for 
the given values of the tail strength \i. Various limiting behaviors can be achieved in eq.(l) 
by varying fx, namely: a) fx » 1 corresponds to P(g) — > S(g — 1), b) fi = 1 to a flattop 
distribution between < g < 1, c) < fx < 1 corresponds to continuum percolation, and 
d) li < to algebraic tails for g — > oo. In addition one can study the limit e) fx — > 
which corresponds to the relevant distribution for dirty semiconductors in the hopping 
regime[13]. 

Some background and notation for transport and percolation exponents is in line: For 
infinite systems near p c above the threshold, the conductivity a behaves as a ~ (p — p c )*, 
where t is the conductivity exponent. In the percolation problem a characteristic length 
£ ~ {p — Pc)~ v of geometrical origin emerges. The correlation length exponent v is defined 
by the previous expression. Below this length scale the system is strongly inhomogeneous 
and has a fractal structure, and thus intensive quantities such as the conductivity scale 
with the length L in a nontrivial manner. Beyond this scale the system becomes effectively 
homogeneous, and intensive quantities take their macroscopic value ( the macroscopic 
average conductance (G) = aL d ~ 2 with a independent of L.) It is an accepted theoretical 
result [7] [8] [9] that when the PDF of conductors corresponds to the continuum percolation 
distribution and the underlying lattice is at p c , one finds a regime where universal lattice 




(1) 
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exponents are found, while values of fi below a certain threshold fi c yield tail strength 
dependent exponents. This behavior is summarized by the expression t = max{to,t(fi)}, 
where to = (d — 2)u + 0o and t(/j,) = (d — 2)u + 

Hierarchical lattices of length L ranging from 2 to 16384 with random microscopic 
conductances obtained from eq.(l), are generated for selected values of fx. The equivalent 
conductance G of each lattice is then calculated exactly. More than 10000 realizations of 
randomness are generated to produce a sample from which a renormalized distribution of 
equivalent conductances can be estimated. Our results are the following: 

1) For lattices above the percolation threshold, we find a rapid convergence of the 
MPD to a sharp Gaussian for all fx ^ 0, as shown in fig. (2). The conductance is then a 
self averaging quantity in the sense that, as the network increases in size the conductance 
approaches a limit value with diminishing fluctuations. For the conductance one finds that 
(G) ~ G typ for large enough L, indicative of a peaked distribution (G typ stands for most 
probable conductance). 

Fig. (3) shows the behavior of the conductance as a function of the length L. Two 
well defined regions are apparent; for small networks the conductance drops sharply, while 
beyond a given disorder length £ D , classical (G ~ L d ~ 2 ) behavior is recovered. The curves 
for different values of /i can be collapsed by scaling L by £d and the conductance G by G^ D 
(the conductance at the disorder length). The result is two universal curves for fi > and 
fi < 0. In the limit fi — > 0, the collapse indicates that £d diverges as a power of the tail 
strength parameter fi i.e. £d = A±\fx\~ z , as depicted in fig. (3). The ratio of the amplitudes 
is A + /A_ = 0.6 and the exponent z = 1.6 ± 0.1. The subscripts on the prefactor indicate 
the approach from either fx > or fi < 0. For lengths larger than £d, the conductivity 
reaches its macroscopic value (given £ < £d). One can then write (in d dimensions) 
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which defines the macroscopic conductivity a = (G^ D / 'Cd~ 2 } ■ One can go further and derive 
the form of (G^ D ) using critical path analysis, valid for £d sufficiently large [13] [16]. This 
is done in the next section. For d e = 3 collapse is achieved by scaling the vertical axis by 
G^ D = (1 — p c ) v+1 / ^ [i v , scaling variable which will be discussed in the next section. In 
figure(4) we show the effective support for conduction on a hierarchical lattice of d e = 2 
conveniently drawn for the purpose of illustration (top of figure). We depict only the bonds 
supporting 99% of the current on the hierarchical network as a function of \i. The disorder 
length £ D derived above is related to the vertical length of the voids. As disorder increases 
(n — > 0) the current is essentially carried by a percolation backbone. 

In summary no matter how broad the MPD (except for n = 0), the limiting behavior 
for P(g) on hierarchical networks of d e = 2, 3 converges to a sharp gaussian. Geometrical 
disorder produced when p c < p < 1 seems not to affect the limit behavior of P(g). A 
disorder length £ D is found for the convergence to the 'Gaussian' basin, which defines the 
length beyond which fluctuations in the conductance decrease around an average value 

(G) = G t yp. 

2) For networks at the threshold, two regimes occur as found by Machta, Guyer and 
Moore (MGM) [17]. For \i > \x c ~ 0.75 the system is self-averaging so that a peaked 
distribution is obtained for large L. In this regime universal lattice exponents are recovered. 
We note that the limit distribution for resistors, although highly peaked, preserves a degree 
of skewness (as measured by the third cumulant) as far as we can determine. This means 
that the limit PDF may not approach a Gaussian behavior as found for lattices above 
the threshold. On the other hand for /j, < ji c the MPD tails are preserved (see fig. (5)) 
and non-universal conductivity exponents are obtained. For the ARC lattice we find a 
disorder length £d by scaling the conductance of the network by L a with a = ln(2/5) / In 3. 
This scale factor corresponds to the scaling of the conductance for the ARC network 
depicted in fig. (2), with no randomness present. Fig. (6) shows evidence of a disorder 
length by collapsing curves for different values of fx. The horizontal axis is scaled by 
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£d = A + (fji — Hc)~ z ' . By collapsing the curves (see fig. (6)) we find fj, c = 0.75 ± 0.05 and 
z' = 3.2 ± 0.1. The value of /j c is close to that reported by MGM on the basis of numerical 
estimates. The exponent z' on the other hand is a new exponent. Presumably this new 
length scale determines a substructure of the ARC lattice supporting most of the current. 
The prefactor A + ~ 500, is large in comparison to the ones obtained above the threshold 
and is indicative of long crossover effects well recognized in the literature[18]. Regarding 
the scaling of the average conductance at the disorder length (G^ D ), we find a very weak 
but clear fi dependence as G^ D ~ B log fx. 

Below /j, c , t; D remains infinite and MDP tails are preserved through the scales as 
reported before[17] [19]. We call this behavior "Levy-like" for its similarities to the case 
of addition of random variables, when the tail strength is such that the first and second 
moments diverge. The additive limit theorem is obviously equivalent to the resistor network 
in one dimension (in terms of resistors). 

II. Disorder Length 

In the last section the existence of a disorder length scale £d for homogenization 
was shown numerically. In the present section we discuss the origin of this length scale, 
for networks above the threshold, starting from the theory of Ambegaokar, Halperin and 
Langer[13]. This theory was originally applied to transport processes that involve quantum 
mechanical tunneling or thermal activation over a barrier, where the barrier distribution 
is itself broad. The relevant distribution for these systems is of the form 

P(g) = ^- with Kg<e\ (3) 

where the range for g is restricted so that the distribution is normalizable. This PDF is so 
broad as A — > oo, that conduction is restricted to a small subset of the whole supporting 
structure. AHL conjectured that the relevant subnetwork would be dominated by the 
bottle-neck resistor which first establishes a conducting path by the following procedure: 
remove all conductors recalling their location of origin, then put them back in order of 
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decreasing conductance. The bond laid at the point of percolation will determine the 
macroscopic conductance. This process is formally expressed by the percolation condition 

poo 

/ P(g)dg = p c , (4) 

•> 9c 

where g c is the governing conductance. This argument is applicable always that A is 
large, so that P(g) in eq.(3) has strong enough tails. In the following we shall establish, 
by plausible arguments, the dependence of the disorder length on the tail strength 
parameter \i. Arguments will pertain mostly to the results above the percolation threshold 
described in section I. 

By definition the correlation length is 

£ = k(p - Pc)- U = k\ v (H9c/9))~ u , (5) 

where we have used eq.(3) and the percolation condition. The prefactor / A^ defines a 
length scale which we identify with the disorder length £d. The argument builds on a 
relation that suggests a use the critical path analysis for more general distributions[16], 
namely 

P(9c)9c = A" 1 « 1, (6) 

which simply states that A should be large. Applying the percolation condition for the 
power law distributions in eq.(l) one arrives at the following dependence for g c , 

g c = pl^ for fi < 0. 

Using the condition in equation (6) yields 

Pc9c = \n\9c = \lAPc = for A» < 0. 

By eq.(5) one can establish that A = {^d/Iq) 1 ^ so for \x < 
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from which 



Analogously for \i > one gets 



(7) 



,(+) 



(1-Pc)" 



From these expressions we can also derive the ratio of the amplitudes (not universal) which 
comes to be 

"(1-Pc 



A+/A. 



Pc 



(8) 



Taking the values for p c and v, which can be computed exactly on these lattices, we 
get A + /A_ = 0.46 in reasonable agreement with the value 0.6 found numerically. This 
simple argument leads to the conclusion that z = v the percolation correlation length 
exponent, and the correspondence \i <-> p. One can also arrive at the same conclusion by 
a simpler derivation as follows: From £ = l (p — p c )~ u and using eqs.(l-4), one obtains 
£ = lo(d^ — 9c )~ u ■ As we are interested in the ji — > limit we expand in powers of n 
and find to first order 



from which the same length scale £o is identified. One can alternatively state the relation 
between the AHL distribution springing from r = roexp(Ae) and the distributions given 
by eq.(l) by choosing e from the appropriate distribution Z)(e)[16]. 

One can readily explain the scaling variable for the conductance in figure (3), by 
invoking a relation first suggested by Ambegaokar, and Kurkijarvi[20] and later discussed 
in general dimensions by Tyc and Halperin[16], 



( g Cd) = Cg c [g c P{g c 



\(d-2)i 



(9) 



where g c is the critical conductance as if A — > oo, P(g) is a broad distribution and v is 
the correlation length exponent for ordinary percolation. From eq.(6) the conductance 
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(G^ D ) depends on A. In d = 2 no A dependence is present and AHL arguments have no 
corrections. As a matter of fact on square lattices it can be shown by duality arguments 
that C = 1 in eq.(9) and AHL arguments are exact. For the general power law distribution 
in eq.(l) one gets, 



valid for \i > 0. Similar formulae (with (1— p c ) — > p c ) hold for \i < 0. In general dimensions 
we expect 



with t = (d — 2)v + l/n the non-universal conductivity exponent, and F\ and F2, universal 
functions. This is precisely the behavior verified by collapsing data in two and three 
dimensions for fx — > 0. When \i is larger, these predictions are expected to deteriorate 
as the volume within £0 gets smaller and the AHL percolation arguments include further 
corrections. 

The previous discussion has only involved networks above the percolation threshold. 
We have no similar understanding for phenomena ocurring at the threshold. Most likely, 
the behavior observed in this range will follow from a renormalization group treatment 
similar to that of Lubensky and Tremblay[8]. This treatment correctly shows (with cor- 
rections due to Machta[9]) that there is a crossover from normal conductivity exponents 
to non-universal ones at \i c for lattices at the threshold. What transpires from the nu- 
merical results is that the percolating lattice with strong disorder for [i> \x c self averages 
after a length scale Above this length scale the system behaves as if no disorder were 
present i.e. the length dependence of the conductivity is purely of geometrical origin. As 
\x approaches /z c , this length diverges giving way to a non-self-averaging situation, which 
explains physically the origin of non-universal exponents. 




for d = 2 
for d = 3 



(10) 




for /i > 
for /I < 



(11) 



9 



III. Discussion 

Usually one says a system is self averaging if one can divide it in subsystems which 
are representative samples of the whole. The size of the subsystems is determined by the 
distance over which correlations persist. When such correlations persist to all scales the 
system is no longer self averaging. In this work we have caracterized how the correlations 
diverge for random resistor networks as a function of the disorder strength as measured 
by the broadness of the MPD. The limit in which correlations diverge depends on whether 
the geometrical support is critical or not. While for lattices above the threshold only expo- 
nentially broad distributions induce diverging correlations, lattices at the threshold show 
infinite correlation at a critical broadness given by ji c . The disorder length ^ D discussed in 
the previous sections is a measure of such correlations. The study of this length scale has 
shown that its critical behavior is analogous to that of ordinary geometrical percolation 
given that p <-> The length scale £d is in fact well known qualitatively in many con- 
texts: In the study of porous media the porosity does not acquire its intensive macroscopic 
value until one measures beyond a given length scale. This effect is expected to influence 
the scaling of other transport related quantities [21] [22] such as the absolute and relative 
permeability of crucial relevance in oil recovery. In the context of ore mining the length 
scale is known as the "nugget effect" and is easily measurable (see ref.[23].) The knowl- 
edge of this length is very important to the scaling theories of these systems, generally 
described in terms of mean field type approaches [24]. One of the most important draw- 
backs of the mean field descriptions is the correct assesment of the fluctuation correlation 
length. On the other hand if by an alternate theory one can access the correct correlation 
length, one could then incorporate such information in the original mean field scheme to 
improve results. One way to do this is to obtain the renormalized distribution beyond the 
crossover length (or merely its mean and variance) and use it as the MPD in the mean field 
self-consistent equation. This procedure has been carried out by MonteCarlo on regular 
lattices [25] with good results. 
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In the following we will briefly discuss other disordered systems that exhibit analogous 
features to RRN with a broad distribution of bond strengths. The simplest disordered 
system with non-trivial properties is maybe a random walker, whose displacement is given 
by a addition of random variables. The limit PDF of a sum of random variables is well 
known to obey the central limit theorem always that the first two moments of the added 
variables exist. When the distribution of the individual variables is broad enough their 
first and second moment diverge the Gaussian basin is not approached. In its place a Levy 
basin arises in which the limit distributions preserve the tails of the underlying random 
variables. In the language of eq.(10) fi c = 2 is the critical value of the tails that separate the 
Gaussian basin from the Levy basin. As implied by our results there is a strong parallelism 
between limit theorems for addition of random variables and RRN for lattices at p c (see 
Section II) where \x c = 0.75, or lattices above p c with [i c = 0. Derrida[26] has established a 
similar analogy to spin glasses. In that model the high temperature phase is characterized 
by a peaked distribution (annealed limit) well described by its average, while the low 
temperature phase (spin glass phase) corresponds to a broad distribution of properties. 
For the Random Energy Model[27](the simplest Spin-Glass) the correspondence n <-> T/T c 
can be established [28], where T c is the critical temperature for the transition. Presumably 
the same will be true for the more complicated Sherrington-Kirkpatrick[29]model. 

Furthermore a similar phenomenon occurs in the context of directed polymers [30] 
where the composition law for random variables is sums of products. Recently Zhang[31] 
found that non-universal exponents occur in the context of directed polymers if one intro- 
duces a broad distribution of random variables. This type of disorder presumably arises, 
in the context of evolving fluid interfaces, due to the power law nature of the breaking 
process that generates a porous medium. This model was also studied by Roux et al[12], 
and the results are very similar to ours; namely the existence of a percolating limit (where 
AHL type arguments apply well) as the broadness of the distribution increases towards 
exponentially wide, and an analog of a 'Levy basin'. Roux et al are also able to find the 
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asymptotic distributions and scaling behavior for the simple limit of the optimal paths. 
It would be interesting to study similar "disorder lengths" and their divergence in these 
closely related systems. 

As a final point we discuss the action of the composition law for random variables 
as a "filter" for broad tails. From the results of the previous sections the stable limit 
PDF can be inferred (by inspection) from the knowledge of the microscopic PDF and 
the manner in which the composition law preserves, reduces or enhances the tail of a 
probability distribution. We will illustrate this point by discussing the simple model of 
Invasion Percolation. In the context of this model, the interplay between composition 
law and PDF yields an explanation for percolating backbone nature of the invading fluid. 
According to the rules of invasion percolation a non-wetting fluid will invade a porous 
medium characterized by a distribution of pore radii Q(r) (or equivalently by a distribution 
of capillary pressures K(P), where r ~ 1/-P), invading one pore at a time and following 
the path of least capillary resistance. For a hierarchical lattice such as the one shown in 
fig.(l) the composition law to produce the equivalent capillary pressure of a cell is 

P(Equivalent) = max\min{P\, P2), min(Ps, Pa)]. (12) 

The Pi stand for the capillary pressures assigned to the bonds of the hierarchical lattice. 
Such a composition law will progressively maximize the tail of any microscopic PDF, re- 
gardless of whether it is broad or not. In fact the former composition law is identical to 
the one of conductances in the limit that MPD is P(g) ~ 1/g [22]. In view of the previous 
conclusion, the AHL criterion holds and the large scale equivalent capillary pressure of the 
porous medium above the percolation threshold is identical to that of the percolating clus- 
ter of smallest capillary pressure. The same argument applies for wetting fluids replacing 
rain — ► max and viceversa in eq.(12). 
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Figure Captions 



Fig. (1 ) Recursive construction of a hierarchical lattice: Each bond at a given generation m 
is converted to a chosen motif, (a) The Berker lattice d e = 2, (b) the ARC lattice of 
d e = log 4/ log 3. 

Fig. (2) Convergence of a broad microscopic distribution towards a Gaussian as L is increased. 

a) Shows the case of the Berker lattice of d e = 2 and ii = 0.4 and b) the Berker lattice 

for d e = 3 and ii = 0.3. The size of the lattices used are: squares L=16, filled circles 

L=256, triangles L=16384 . The inset shows a fit to as Gaussian distribution. 
Fig. (3) The figure shows the conductance versus the scaling variable L/£ D for a lattice of 

d e = 2, above the threshold. A crossover length £d and a characteristic conductance 

G^ D are identified by collapse of different \i values. The length £ D diverges when 

/I -> as A±\fj,\- 1 - 6±0 - 1 , where A + /A_ = 0.6. 
Fig. (4) Lattice of bonds carrying 99% of the current as a function of fx. The hierarchical 

lattice generator is conveniently drawn as shown in the figure. Only the lengths along 

the vertical axis in the figure are meaningful 
Fig. (5) Limit distribution for a percolating geometry (ARC lattice) and ii = 0.5 < ii c . The 

resulting distribution shifts towards low G values, while power law tails are preserved. 

Log-log scales are used so power law tails are manifest. 
Fig. (6) Conductance curves versus system length for the ARC lattice with fi > fi c . The 

collapse is achieved using £ D = A + ((i — fi c )~ z \ where z' = 3.2 ±0.1 A + ~ 500. On 

the vertical axis a = ln(2/5)/ln3. 
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